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1 INTRODUCTION 



In our recent works lVranies fc Poedtj l|2009al lbll3ldh the drift 
wave theory has been applied to the problem of the heat- 
ing of the solar corona. The drift wave is driven by gradi- 
ents of the background plasma parameters (e.g., gradients 
of the plasma density, temperature, and magnetic field) that 
may act either separately or all together. In the literature, 
the drift wave is frequently called the universally unstable 
mode because it is growing both within the (two- or multi- 
component) fluid theory and the kinetic theory. In the for- 
mer (fluid) description, it is either a coUisional instability 
(discussed iniVranies & Poedts (2006) in the context of the 
solar atmosphere) that appears as a combined effect of the 
electron collisions, the ion inertia, and the density gradient, 
or a reactive one (|Vranies fc Poedtsll2009dl ) in the presence 
of all three gradients mentioned above where now the ions 
play a major role. The latter (kinetic) theory, on the other 
hand, can also describe these two instabilities. However, the 
kinetic description yields an additional (purely kinetic) in- 
stability that is absent within the fluid description. In appli- 
cation to the solar cor ona, this additi onal kinetic instability 
is studied in iVranies fc PoedtsI l|2009al la). This purely kinetic 
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ABSTRACT 

In the solar corona, several mechanisms of the drift wave instability can make the mode 
growing up to amplitudes at which particle acceleration and stochastic heating by the 
drift wave take place. The stochastic heating, well known from laboratory plasma 
physics where it has been confirmed in numerous experiments, has been completely 
ignored in past studies of coronal heating. However, in the present study and in our 
very recent works it has been shown that the inhomogeneous coronal plasma is, in fact, 
a perfect environment for fast growing drift waves. As a matter of fact, the large growth 
rates are typically of the same order as the wave frequency. The consequent heating 
rates may exceed the required values for a sustained coronal heating by several orders 
of magnitude. Some aspects of these phenomena are investigated here. In particular 
the analysis of the particle dynamics within the growing wave is compared with the 
corresponding fluid analysis. While both of them predict the stochastic heating, the 
threshold for the heating obtained from the single particle analysis is higher. The 
explanation for this effect is given. 

Key words: Sun: corona, oscillations. 



growth of the drift wave is entirely due to electron dynamics 
and it develops provided the wave frequencies are below the 
electron diamagnetic drift frequency. 

Regardless which of the mentioned instabilities is in 
action, the result is always a growing wave amplitude. 
When some critical value for the wave potential is ac h ieved , 
stochastic heating sets in (JMcChesnev et all 1 19871 . Il99ll : 
ISanders et al.l [l998J). By nature this heating is strongly 
anisotropic (acting mainly in the direction perpendicular 
to the magnetic field vector). Moreover, it is mass depen- 
dent (heating mainly ions, with better heating of the heavier 
ions) , it is due to the electrostatic features of the drift wave, 
and it is a low frequency process (typically the drift wave 
frequency is much below the ion gyro-frequency). 

The drift wave can easily become electromagnetic in 
case of a relatively high plasma-/? (exceeding the electron 
to ion mass ratio). In that limit, it is in fact coupled to 
the Alfven wave. Yet, even in this case the heating de- 
velops almost without change and primarily due to the 
electric field of the drift- Alfven mode. In the past, this 
has been studied both analytically and numerically, and it 
was even experimentally verified in laboratory (tokamak) 
plasmas llMcChesnev et aLlll987l . Il99ll : ISanders et aLlll998l : 
I White et al.ll2002l ). For the coronal environment the electro- 
magnetic drift-Alfven ca se has been discussed in detail in 
IVranies fc PoedtsI l|2010l ). 
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For the purpose of the present study, in order to quanti- 
tatively describe some details of the heating and acceleration 
of plasma particles, any of the mentioned instabilities of the 
drift wave will do. Hence, we shall use the kinet i c, dens ity 
gradient driven instability from lVranies fc Poedtd (|2009al lbr). 



THE INSTABILITY OF THE DENSITY 
GRADIENT DRIVEN DRIFT WAVE 



From standard textbooks (jlchimarul [l980l : IWeilandll2000t ). 
the frequency and the growth rate of the kinetic drift wave 
are given by 
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Here, the geometry is chosen as follows. The constant mag- 
netic field is given by Bq = BqEz, the equilibrium density 
gradient in the electrically neutral plasma with nm = rieo = 
no is given by Vxno = —exn'o = —Exdno/dx, and the per- 
turbations are assumed to be of the form f{x) exp(— jojf + 
ikyy + ikzz), with \kz\ <^ \ky\ and an x-dependent wave am- 
plitude /. The local approximation will be used. Hence, it 
will be assumed that |d/da;| <^ \ky\. 

In the derivations of Eqs. ((TJ and[2|, the condition of 
a strongly magnetized plasma is used, \ij\ ^ Qi, together 
with the smallness of the acoustic response of ions in the 
direction parallel to the magnetic field vector, implying that 
Ur/kz S> Cs, where Cs is the ion sound speed. The notation 
in Eqs. (O and (0) is standard (|Weilandll2o"ool '). i.e.. 



= -fc. 



2 / 



-{Ti/Te)uj*e, VT.j=KTj/mj, 



Ao(bO = Io{bi)eyip{~bi), bi = kypi, \di = Vri/uJ^i, 

where 7o is the modified Bessel function of the first kind and 
of the order 0. 

From Eq. ((2)| it is seen that within the present insta- 
bility, the sign of 7 can be changed only by the electron 
term if \ujr\ < \iiJ,e\. This is, however, only a necessary con- 
dition. A sufficient condition is obtained when the contribu- 
tion of the rest of the expression (the ion damping terms) is 
taken into account. As d emonstrated by various graphs in 
IVranies fc PoedtsI (|2009al la) , in the case of the solar corona it 
is more difficult to find a regime where the mode is damped 
than a regime where it is growing. In these studies it is 
shown that the mode grows more strongly for shorter per- 
pendicular wave lengths Xy (meter or sub-meter size) and 
for longer parallel wavelengths Xz that may exceed the per- 
pendicular ones for 4 to 5 orders of magnitude. Note that 
a similar ratio of the two components Xj holds for tokamak 
plasmas as well. The same conclusions can also be drawn 
in the ca se of the reactive (the so-c alled jyi-instability) dis- 
cussed in lVranies fc PoedtsI l|2009dr ). 

To provide some details of the particle acceleration and 
heating by the growing drift wave in the context of the solar 



corona, we need specific data. Hence, as an examp l e we he re 
use a set of parameters from IVranies fc PoedtsI (|2009al lbl) 
that are known to yield the instability: Bo = 10"'^ T, 
no = lO^'^ m'^ L„ = [{dno/dx)/no]~'^ = s ■ 100 m, 
Xy — 0.5 m, and Xz — s ■ 10'* m. These are just representa- 
ti ve parameters with no part icular significance. As discussed 
in lVranies fc PoedtsI l|2009br ). the instability will in fact de- 
velop even when the density is varied by several orders of 
magnitude around the given value, and the same holds for 
the magnetic field and the tempe r ature. The parameter s is 
introduced in IVranies fc PoedtsI (|2009al lbr) for convenience 
only, because it was realized that the ratio y/tUr remains ex- 
actly the same in case the ratio Xz/L„ is fixed. Here, L„ is 
the characteristic scale-length for the inhomogeneous equi- 
librium density. This further implies that we may go to very 
different scales, e.g., by taking s in the range 0.1 — 10*. As 
long as Xy is kept constant, the ratio ^/ujr will remain the 
same, although the actual values for these two quantities 
will certainly change. Hence, for the given parameters, from 
Eqs. ([ij and ((2)1 one finds j/uj,. = 0.26. In the case s = 1, 
this in fact implies Ur — 254 Hz and 7 = 66 Hz. While in 
the case s — 10"^, this would give a mode with LUr — 0.254 Hz 
and 7 — 0.66 Hz. Going to such small values of Ln (tens of 
meters) in fact makes sense because the perpendicular diffu- 
sion, that naturally develops in such a plasma with density 
gradients, is indeed very small, with the diffusion coefficient 
D± of the order of 0.01 m^/s and the corresponding diffu- 
sion velocity D±V ±n/n of just a few millimeters per second 
(|Vranies fc PoedtauOOSi ). Short values of L„ in the same 
time imply very high frequency drift waves and fast grow- 
ing instabilities, so that the expected heating will develop 
at time scales that are orders of magnitude shorter than the 
plasma diffusion. Note that the parameters Xz and L„ can 
of course also be varied independently of each other. 



3 PHYSICAL MECHANISM OF HEATING BY 
THE DRIFT WAVE 

3.1 Heating within the collective interaction 
frame 

The physics of the mentioned stochastic heating is described 
in detail in vario us sources. He re, we give some short sketches 
of it by following iBellanI (|2006l ). and using the standard drift 
wave theory. Within the two-fluid theory, the perpendicular 
velocity of the ion unit volume is give n by the recurrent 
formula (jVranies fc Poedtsll2006l.l2009bD 

1 ^ ^ , «Tj - V±ni V± • TTi 

Vii_ = -p-ez X V±<^ -I- -pr-Bz X h e^ X — 

-DO S2i Hi nuniili 

Id^^ d d ^ „ 

The four terms here correspond to, respectively, the E x B- 
drift that is usually the leading order one, the diamagnetic 
drift w*;, the stress tensor drift u,ri, and the polarization 
drift Vpi. The velocity can be calculated up to small terms 
of any order using the drift approximation \d/dt\ -C ^i- 
In some situations, in plasmas with relatively cold ions, as 
a first approximation one may assume the E x B-drift as 
the leading order one, and set it into the velocity Vi± in the 
polarization drift. The diamagnetic drift does not contribute 
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to the ion flux V ■ (n««i) = 0, and also its contribution to 
the convective derivative of the polarizati on drift is canceled 
by the part of the stress tensor flux nvni ijVranies fc Poedtj 
l2009el ). Therefore, we shall focus now on the E x B-drift, 
and the polarization drift Vpi = {dE ± / dt) / {BoCli) , where 
E±^ = —eyky(j>i. These two together describe an elliptical 
motion in the presence of the wave, described by: 



ky(j)l 



\{kyy + kzZ — curt), 



(4) 
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In deriving Eq. (O the wave amplitude is taken as time 
dependent due to the wave instability, cj>\ = (f> e-icpi^t) , 
cf> — const., and the same holds for the coordinates y{t),z{t) 
appearing in the wave phase ip{t). That is the reason for 
the appearance of the terms with 7 and dz/dt. Neglecting 
these two terms, and taking a constant potential, Eq. (O be- 
comes identical to the corresponding expression from lBellanI 
(1200611 . 

Crucial for the heating is the term 
[kyCJ>i{t)/{BQQ,i)]cos(p in the denominator. Looking 
back into the derivations, it turns out that it appears from 
the derivative d{simf)/dt — ^ dy/dt = Vy = Vpi. Keeping it 
here has sense only in case of a relatively large displacement 
in the direction of the perpendicular wave vector. From 
Eq. ^ it is seen that for 



a{t) = fc^0i(t)/(Bo«,:) = klpte<f>i{t)/{t^T,) ^ 1, 



(6) 



the denominator may vanish and, consequently, the polar- 
ization drift should go to inflnity. Note however, that Eq. ((5]) 
is obtained approximately starting from the guiding center 
approximation that becomes violated for such a large po- 
larization drift, and the analytical expression for the polar- 
ization drift in this limit becomes inaccurate. Yet, Eq. Q 
describes the general trend for the polarization drift; a fast 
and signiflcant acceleration and heating of the plasma par- 
ticles will definite l y take place , as experimentally ver ified in 
iMcChesnev et al.l l|l987l . ll99ll ): [San"ders et all l|l998ll . 



3.2 Individual particle dynamics 

As discussed in iBellanI (|2006l ) , a direct numerical integra- 
tion for plasma particles is required for a full validation of 
the above described heating mechanism for the solar corona. 
The numerical simulation has to take into account the ac- 
tual values of the electric field in the position of a specific 
gyrating plasma particle. For that purpose, taking again the 
drift-wave electric field with a time-varying amplitude, the 
appropriate set of equations for a single particle reads 



y"{t) + y{t) -a(i) sin 
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y{t) + -^zit) - ht 

Ky 



ht 



= 0, 



0, 



(7) 

(8) 
(9) 



0.015 



0.000 



-0.015 




-120 



Figure 1. Particle positions in the x, y-plane within a wave pe- 
riod, after it starts from the point A with (x, y, z) = (0, 0, 0), for 
a small wave amplitude </)i = 0.86 V. 
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Figure 2. Normalized particle positions in the x, y-plane after 
it starts from the point A with {x, y, z) = (0, 0, 0), for a growing 
electric field potential 4>i{t) = (j)exp{jt/Qi), and <p = 0.86 V. 



Here, the prime denotes a derivative in time, and the 
spatial and time coordinates are normalized as x,y,z -^ 
kyX, kyy, kyZ, t — >■ flit, b = iUr/^i- Equation (|9} is obtained 
after one integration where an integration constant appears 
that has been neglected. Note that Eq. (O can be reduced 
to the Mathieu equation with a source term. Namely, ne- 
glecting the z-terms and after introducing r = bt/2, for the 
small remaining argument ip it can be written in the generic 
form 



d y/dr + [a — 2q{t) cos2r]j/(r) = c{t) sin2r. 



(10) 



where a = 4/6^, q(r) = 2a(T)/6^, and c{t) = -4a(T)/6^. 
The solutions of this equation without the right-hand side 
are the Mathieu functions M(a, q, r). Equation (|lUp has un- 
stable solutions and it is discussed in detail in ?. I White et al.l 
(2002). 

The set of ordinary equations ((Zll-® is solved numer- 
ically. We do this first for a very small and constant wave 
amplitude in order to demonstrate some basic features of 
the particle motion due to the wave electric field. The result 
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Figure 3. Normalized velocities Vx,Vz for the same parameters 
as in Fie. 2. 
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Figure 4. Change in time of the normahzed kinetic energy of a 
particle with a unit mass E = (ti^ + u^ + v'^)/2. 



is shown in Fig. 1, where we give the projection of the pro- 
ton positions in the perpendicular plane, x{t),y(t), in the 
wave field with the constant amplitude e<j)\/{KTi) = 0.01, 
that corresponds to <^i = 0.86 V. Here and further in the 
text we take \y = 0.5 m, A^ = 2 • 10* m, which consequently 
from Eq. ([l| yields the wave frequency ujr ~ 254 Hz. The 
normalization is the same as above, viz. {x,y) = {kyX,kyy). 
This all yields a = 0.014, hence the stochastic heating is ab- 
sent, and the particle trajectory is a very elongated ellipse, 
with the motion almost completely in the a;-direction due 
to the E X B-drift described above (observe the difference 
in I, y-scales). After a wave period the particle returns to 
the starting position A with the coordinates {x,y) — (0, 0). 
The ratio of the two axes of the ellipse sijei = 0.00013 de- 
scribes the fact that the particle motion is mainly due to 
the E X _B-drift in the a:-direction, and that the polarization 
drift in the direction of the wave-vector (the y-direction) is 
negligible. 

Performing the same plot for electrons yields the same 
€2 as expected (the E x _B-drift is independent on mass), 
while ei becomes reduced by the factor m^/rrii (the polar- 
ization drift separates charges and it is mass dependent). 



On the other hand, taking a time-varying potential 
<^i(t) due to the above demonstrated drift wave instability, 
and restricting the time interval to relatively short values 
(below the onset time of the stochastic heating), yields a 
spiral trajectory in the x, y-plane. One example of this is 
presented in Fig. 2. Here, the starting value for the poten- 
tial is the same as above e4>i{Q) / [KTi) — 0.01, and the same 
holds for the particle position. The other parameters are the 
same as in Fig. 1. With the same normalization the potential 
is given by (/f>i(i) — 4> expi^'yt / Q,i) , 7 = 0.26cJr, tJr = 254 Hz, 
<j) — <?f>i(0). The maximum time in physical units here is 
0.04 s, and in this moment the potential has reached the 
value of 17 V. Comparing with Fig. 1, one observes that the 
displacement in the y-direction due to the polarization drift 
is now increased by a factor 20. 

It is interesting to compare the leading perpendicular 
_Exi3-drift velocity Vx, and the parallel velocity v^ (along the 
magnetic field vector). The parametric plot in Fig. 3 shows 
that for these potential amplitudes the parallel velocity Vz 
remains about one order of magnitude smaller. 

In order to see what happens for larger amplitudes of 
the electrostatic drift wave when the stochastic heating is 
supposed to be in action, Eqs. ©-([Ql) are solved by allowing 
a slightly larger time range. One of the results is shown in 
Fig. 4 for the kinetic energy of a particle with unit mass 
E{t)/m = [vo,{t)'^ + Vy(t)'^ + Vz(€f\l1, with normalized ve- 
locities as above. The other parameters are the same as in 
the previous text and figures. Here, the stochastic heating 
takes place after around 0.078 s, in the moment when the 
growing wave amplitude reaches the value of around 150 V. 
Note that this is by about a factor 2.5 larger than the value 
obtained from the condition ((6| which, in fact, follows from 
an approximative procedure as explained earlier. Thus, a 
higher necessary threshold for the stochastic heating is ex- 
pected. However, a plasma can support multiple waves in 
the same time. The drift wave spectrum described by Eqs. 
([T][2| in realistic situations imply the presence of more waves 
rather than a single one. The analysis presented in a recent 
study (jSheng et al.ll2009l ) shows that in such cases the insta- 
bility threshold can be considerably reduced. Therefore, the 
ion heating by the mechanism which we discuss here will be 
even more efficient. 

The corresponding plot of the displacement in the di- 
rection of the wave j/(t), which is in fact also equal to the 
velocity v^if) [see Eq. ([9])] is shown in Fig. 5. The plot of 
the velocity Vy (t) (that is not given here) reveals the ampli- 
tude that is completely negligible until Q,it becomes close to 
8 • 10*; after that it is very stochastic and with the ampli- 
tude equal to that of the velocity Vx{t). Thus we are in the 
range of parameters in which the expansion procedure used 
previously is not valid any longer because the polarization 
and E X B drifts are of the same order, and the stochastic 
heating is fully in action. Note that, although the motion 
is stochastic/chaotic, it is in fact de terministic and not ran- 
dom as commented in lBellanI (|200y). Observe that the total 
physical time interval for Figs. 4 and 5 is only about 3-4 
wave-periods. 

The plot of the parallel velocity Uz(t), presented in 
Fig. 6, shows that along the magnetic field the particle dy- 
namics remains non-stochastic, as should be expected in 
view of the analysis presented in the text above. The ve- 
locity magnitude is below 1, thus still smaller than the per- 
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Figure 5. Normalized displacement in the direction of the po- 
larization drift j/{t), and the normalized perpendicular velocity 
component j;(4) = v^it). 
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Figure 6. Normalized perturbed velocity in the direction parallel 
to the magnetic field. 



pendicular velocity (see Fig. 5) which is the consequence of 

the stochastic heating. 

The theory pres ented in lMcChesnev et all (|l987l . ll99il ): 
ISanders et al.l (|l998l ) predicts the maximum stochastic ve- 



locity Vmax — [kyPie(j)/iKTi) + 1.9]Qi/ky, which further 
yields the effective achieved stochastic temperature Tmax = 
mi)^^3./(3K). For the parameters used above this yields 
Tmax — 2 • lO'' K, i.e., the increase of the temperature for 
one million Kelvin within the growth time Tg ~ 0.1 s. The 
maximum energy released by this mechanism is "Smax = 
norriiVmax = 0.04 J/m~^. The energy rele ase rate, Smax/Tg, 
exceeds for several orders of magnitude IjVranies fc Poedtj 
l2009ai ibi') the value presently accepted by researchers as re- 
quired for a sustained coronal heating. 



4 SUMMARY 

The present study provides some additional details of our 
recently proposed model for the coronal heating, which im- 
plies the heating by electrostatic drift waves. The heat- 
ing by waves, which presently represents one of two ma- 



jor heating scenarios, has been studied in many works from 
the early d ays o f the coronal hea t ing in vestigation, e.g., 
by lAlfv6nl (Il947l') and [ Piddington fl956'), but also more 
recently in IPekiinlii et al. ( 2001 ); Suzuki ( 2004 ). More de- 
tails about t he existing heating models are available in var i- 
ous sources (JNarain fc Umschneideij[l99(]| : lKlimtchukll2006l ) , 
some of the m have also been d i scussed in more detail in our 
recent work lVranies fc PoedtsI l|2009bl ). 

In this work, the trajectories of individual plasma par- 
ticles have been investigated within the framework of the 
drift wave instability and the associated stochastic heating, 
with the aim to describe some fine details of the particle 
dynamics during the heating process. The calculation of the 
actual growth rates of the mode is not repeated here because 
it has been given in detail in our recent studies. From the 
present study and from our recent works it follows that the 
general trend of the ion behavior in the wave field with a 
growing amplitude can be described analytically from the 
fluid point of view, as discussed here in Sec. 3.1. However, a 
more detailed picture is obtained by following a single par- 
ticle which is moving in the wave field. For example, this is 
seen from Sec. 3.2 in the case of the critical threshold for 
the onset of the heating. We observe that the general theory 
of the drift wave instability presented in some early studies 
(JTerashimal 1 19671 ) . and not necessarily related to the heat- 
ing, also includes the particle aspect approach, and as such, 
is in perfect agreement with the fluid and kinetic modeling. 

It is also interesting to note that, although the individ- 
ual particle dynamics becomes stochastic for relatively large 
wave amplitudes, the overall picture of the macroscopic mo- 
tion of the plasma, caused by the wave, actually remains well 
described by the fluid d escription. This has been directly ex- 
perimentally verified bv lBailev et al.l (| 19931 ). The parameter 
a [c.f. Eq. ((B))] in their experiment was around 2, and thus 
the stochastic heating was in action, yet the observed elec- 
trostatic potential and density profiles nicely matched the 
theoretical contour plots corresponding to a wave described 
within the fluid theory and drift approximation. This is in 
agreement with the statement given in Sec. 3.2 about the 
nature of the heating; the increased effective temperature 
Tmax follows froiu the stochastic velocity Vmax which in fact 
is not random, yet in practical terms it is indistinguishable 
from standard coUisional thermalization and the broadening 
of the distribution function (|Bellan 2006i ). 

The model used here implies a maxwellian distribu- 
tion function for plasma species, that is the basis for 
both the kinetic and the fluid description of the instabil- 
ity. However, t he stochastic heating with all its features 
(|Vranies fc Poe dts 2009b d) in fact implies that such a start- 
ing distribution function will necessarily evolve. The perpen- 
dicular ion temperature will grow and eventually become 
larger than the parallel one. Hence, the particle distribu- 
tion may pass through several stages, going through the 
Oort's distribution with the temperature anisotropy and to 
the more general loss-cone distribution (jBaronia fc Tiwaril 
1999), or to something else. This opens an interesting pos- 
sibility for eventual numerical simulations because such fine 
details are beyond the scope of an analytical study. 
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